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Modeling color-dependent galaxy clustering in cosmological 
simulations 
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ABSTRACT 

We extend the subhalo abundance matching method to assign galaxy color to subhalos. We 
separate a luminosity-binned subhalo sample into two groups by a secondary subhalo prop- 
erty which is presumed to be correlated with galaxy color. The two subsamples then represent 
red and blue galaxy populations. We explore two models for the secondary propertty; sub- 
halo assembly time and local dark matter density around each subhalo. The model predictions 
for the galaxy two-point correlation functions are compared with the recent results from the 
Sloan Digital Sky Survey. We show that the observed color dependence of galaxy clustering 
can be reproduced well by our method applied to cosmological TV-body simulations with- 
out baryonic components. We then compare the model predictions for the color-dependent 
galaxy-mass cross correlation functions with the results from gravitational lensing observa- 
tions. The comparison allows us to distinguish the models, and also to discuss what subhalo 
property should be used to assign color to subhalos accurately. We show that the extended 
abundance matching method using local dark matter density as a color proxy provides an 
accurate description of the galaxy populations in the local universe. 

Key words: galaxies: evolution; galaxies: haloes; galaxies: statistics; cosmology: large-scale 
structure of Universe; gravitational lensing: weak 



1 INTRODUCTION 

Describing galaxy formation and evolution within the context of 
the standard cosmology is one of the most important goals in as- 
tronomy and cosmology. Large data from recent gala xy redsfhit 
surve ys, such as the Sloan Digital Sky Survey (SPSS: lYork et al] 
l200Cb . have enabled us to measure precisely various statistics 
for subsamples of the galaxies classified by, for example, lumi- 
nosity, stellar mas s, morphology and so on (e.g., iLi et al. I l2006l : 
IZehavi et alfeOllh . The upcoming deep and wide-area galaxy sur- 
veys, e.g., Subaru Hyper Suprime-Cam Survey and Dark Energy 
SurvejQ will provide even larger data of distant galaxies. Clearly, 
accurate theoretical models are needed to interpret observational 
results such as those on the relationship between galaxies and dark 
matter (DM) halos. 

So far, several useful phenomenological methods to 
link halos with galaxies have been developed. One is the 
so-called halo occupation distribution (HOD) modeling 



* E-mail: shogo.masaki@nagoya-u.jp 

f JSPS Fellow 

1 http://subarutelescope.org/Projects/HSC/ also http://sumire.ipmu.jp/en/ 

2 http://www.darkenergysurvey.org/ 



(see e. g.. ISeliakl l200d: ICoorav & Shethl 120021: iTinker et al l 



20051: iBlake. Collister. & Lahavl 120081 ; IZehavi et all |201ll 



Leauthaud et all 20121) . The method parametrizes HOD, i.e., 
occupation number of galaxies per halo, of central and satellite 
galaxies separately, as a function of halo mass. The functional form 
of HOD is inspired by hydrod ynamical cosmolog ical simulations 



ot HUJJ is inspired by nydrodynamical cosmological simulations 
or by semi-analytic models dZheng et al.1 12005 ). By using the 
halo mass function as well as the bias factor and the DM density 
profile obtained from cosmological TV-body simulations, one can 
accurately calculate the galaxy two-point correlation functions. 



There are free parameters in the HOD formalism which 
are constrained by matching to the observed clustering and 
to th e number density o f galaxies 1 Blake. Collister. & Lahavl 
120081 : IZehavi et al.l 1201 it ICoupon et alj l2012h . More sophisti- 
cated treatments of HOD include the con d itional lumino sity func- 
tions jvan den Bosch. Yang. & Mcl 120031 : ICooravl l2006h . An at- 
tractive feature of HOD is that it can be extended to model 



and calcu late other statistical quantities than galaxy cluster- 
ing (e.g., De Bernardis & Coorav 2012 : Leauthaud et al.1 120121 



iMasaki & Yoshidall2012l : Ivan den Bosch et alj|2012l) . However, it 
is known that one often encounters difficulty in fitting HOD 
model predictions to the observed clustering and the number 
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density simultane ously dOuadri et al.ll2008l; iMatsuoka et ai]|201ll ; 
IWake et aljioTTh . 

The second method is called the subhalo abundance 



20041; Conrov. Wechsler. & Kravtsov 


20061; iGuo et al. 


20 id; 


Moster et alj 1201 & 1 Wetzel & White 


l20ld: iNeistein et al.l 


20111; iTruiillo-Gomez et alj 1201 it iReddick et al.l 


20121: 


Rodn'guez-Puebla, Drorv. & Avila-Reesd 2012b. SHAM as- 



sumes a tight relation between galaxy luminosity (or stellar mass) 
and subhalo circular velocity (or mass). The threshold circular 
velocity is set such that the number density of subhalos and that of 
luminosity-selected galaxies are matched as 



"-subhalo (> Vdi) = ?l ga l axy (> L) . 



(1) 



The predicted clustering generally agrees with observations, repro- 
ducing the observed luminosity dependence of galaxy clustering. It 
should be emphasized that SHAM does not need fitting parameters 
while HOD modeling does. An additional advantage of SHAM is 
that it can utilize iV-body simulation outputs directly. Therefore, 
SHAM can incorporate the non-linear evolution of galaxy clus- 
tering accurately while having a simple and direct galaxy-subhalo 
connection. 

Both the HOD model and SHAM are used to calculate the 
galaxy two-point correlation or the galaxy-mass cross correlation 
as a function of a galaxy property, e.g., luminosity and stellar mass. 
It would be ideal if multiple observable properties are assigned si- 
multaneously to (sub)halos. Galaxy color is one of the most fun- 
damental properties. It is thought to indicate the galaxy's age and 
the level of star formation activity. Also it is well-known that there 
are two apparent sequences in the color-magnitude diagram, which 
represent red and blue galaxies. The two populations differ in many 
aspects including their spatial clustering, reflecting probably their 
formation histories. So far, detailed treatments for color-dependent 
clusterings or lensing m easurements have been developed, but only 
within the HOD model I van den Bosch, Yan g. & Mo] l2003l;ICooravl 



2006; 



2009 



Mandelbaum et alj2006l;lTinker et al.l2008l;lRoss & Bninnen 
Simon et alj2009l ; ISkibba & Shethl2009l ; IZehavi et al.l201 lb . 



In the present paper, we develop a method to assign galaxy 
color to subhalos by extending SHAM. We first apply SHAM to a 
subhalo catalog to obtain luminosity-selected subhalo samples. We 
then divide subhalos in a luminosity bin into two groups by order- 
ing subhalos using a secondary quantity. The two groups are then 
meant to represent red and blue galaxy samples. Similarly to the 
original SHAM, the abundance ratio of the divided two groups is 
matched to the observed red/blue ratio. Clearly, we need to choose 
an appropriate subhalo property which is presumably correlated 
with galaxy color. We propose two models for the secondary prop- 
erty. One is motivated by the so-called assembly bias whereas the 
other incorporates environment effects. 

Assembly bias is a property-dependent bias for halos in 
a fixed mass bin. A variety of properties have been consid- 
ered so far, such as assembly time, spin, concentration and so 
on (see e.g.. iGao. Springel, & Whitell2005l; IWechsler et al.l | 2006| ; 
Croton. Gap. & White 2007; lGao & Whitel2007l ; lReed et alfeOOTI ; 



Lacerna & P adilla 20111 120121) . Among these, we choose subahlo 



age as a proxy of galaxy color. It is thought that there is some en- 
vironm ental effect through which a galaxy's morphology and color 
evolve. iGerke et alj d2012l) recently developed a model to assign 
color to subhalos (see also iTasitsiomi et alj|2004l) . In their model, 
color is assigned to subhalos by using local galaxy d ensity along 
with the empirical color-galaxy density relati on (e.g., iHogg et aU 
l2004l ; lBalogh et alj| 20041 ; lBamford et aUl2009l) . Our second model 



is motivated by their approach, but we do not use such an empirical 
relation. 

Finally, we consider lensing measurements in addition to 
galaxy clusterings in order to distinguish the two models. Lensing 
observations will provide an independent information which possi- 
bly allows us to derive an accurate galaxy-halo connection. In par- 
ticular, we utilize the observed e arly/late-type dependent lensing 
profile dMandelbaum et al]|2006l) assuming early /late-type galax- 
ies approximately correspond to red/blue galaxies. We show, for 
the first time, that an extended SHAM can reproduce the observed 
lensing profiles as well as the color-dependent galaxy clusterings 
in the local universe. 

The rest of the present paper is organized as follows. In Sec- 
tion 2, we describe our cosmological iV-body simulations and the 
basics of SHAM. In Section 3, we explain our method of galaxy 
color assignment. Then in Section 4, we compare the model predic- 
tions for color dependent galaxy clustering and galaxy-mass cross 
correlations with the observational results. We give summary and 
discussion in Section 5. 



2 SETUP 

2.1 Cosmological iV-body simulations 

We use the massively parallelized co de Gadget-2 
dSpringel. Yoshida. & Whitdl200ld ; ls"pringell2005l) in its Tree-PM 
mode to run two cosmological iV-body simulations. For each 
run, the assumed A-cold-DM cosmology is consistent with the 
WMAP 7-year results, with Q m = 0.272, 6 = 0.0441, fi A = 
0.728, #0 = WOh = 70.2 km s' 1 Mpc" 1 , a 8 = 0-80 7 and 
n 3 — 0.961 in the standard notation dKomatsuetal]|201lh . The 
two runs, one with 200 /i -1 Mpc on a side and the other with 
300 ft -1 Mpc box size, both employ 1024 3 DM particles, and will 
be referred to as L200 and L300, respectively. The resulting mass 
of aparticleis 5.6 x 10 8 h~ 1 M (s (1.9 x 10 9 fe~ 1 M P) ) for the L200 
(L300 ) run. We use the code CAMB dlxwis, Challinor & Lasenbvl 
l200fj|) to obtain the matter power spectrum for the initial condition. 
We set the initial redshift to be 49 (65) for the L200 (L300) run. 
The gravitational softening length is set to be 5 /i _1 kpc and 
8.8 h^kpc for the L200 and the L300 run, respectively. The L300 
run has lower mass and spatial resolution than the L200 run but 
provides higher statistical precision of correlation functions for 
massive objects and at the large scale. 



2.2 Identification of halos, tracking assembly history and 
construction of subhalo catalog 

We identify the distinct halos, i.e., halos that do not lie in more mas- 
sive halos, using the friends-of-friends (FoF) algorithm with linking 
parameter of 0.2 in units of the mean interparticle separation. For 
the identification of subhalos, i.e., dense self-gravitating clumps 
that reside i n a distinct halo, we u tilize the SubFind algorithm de- 
veloped by ISpringel et all d2001al) . The algorithm decomposes a 
distinct halo into a central subhalo, i.e., the so-called smooth com- 
ponent which contains the majority of mass, and satellite subhalos. 
We store the subhalos with more than 20 particles. 

To implement SHAM, we construct the subhalo catalog with 



value of particle circular velocity V c ir(Y) = GA1(< r)/r 
at the last epoch when the satellite subhalo was a central one 
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dConrov. Wechsler. & Kravtsov 2006|). We follow the recipe de- 
scribed in detail by lAUgood d2005n " to construct the mass assem- 
bly history of the most massive progenitor of subhalos. We use 50 
simulation snapshots taken evenly in ln(l + z) from z = 10 to 
2 = for the two simulations. For the central subhalos, we tab- 
ulate the maximum circular velocity at the redshift when they are 
observed, KUax • The term Vdi in Equation |QJ is replaced with 
Knax f° r the central subhalos and V^ax f° r the satellite subhalos. 
Note that the circular velocity defined in this way is a direct proxy 
of gravitational potential and is less sensitive to subhalo identifica- 
tion algorithms. 



2.3 Application of SHAM 

Following lZehavi et alj J201 lh . we consider galaxies in three mag- 
nitude bins, -22 < M r - 51og 10 h < -21, - 21 < M r - 
51og 10 /i < — 20and— 20 < M r — 51og 10 /i < — 19, and hereafter 
denote the absolute magnitude just as M r without the — 51og 10 h 
term. We apply SHAM to the subhalo catalogs constructed in the 
L200 and L300 runs to create the magnitude-binned subhalo sam- 
ples. We take the bracketing threshold samples for a bin and then 
use the difference of them as a binned sample. The co moving num- 
ber de nsity for each threshold sample can be found in lZehavi et all 
42011b . 

We measure the projected correlation function w p as a func- 
tion of the projected distance r p as 



w p (r P ) = 2 / djv £( r = yWrf). 



(2) 



We take the same value of 7r ma x used by IZehavi et alj 4201 lh for 
each bin. For the fainter and the intermediate samples we use the 
L200 subhalo catalog at z = 0. Only for the brighter sample, we 
use the catalog from the L300 run at z = 0.1. The chosen red- 
shift is very close to the mean redshift of eac h magnitude binned 
galaxy sample. The results are compared wifh lZehavi etal.ld201lh 
in Figure [T] The top panel shows fairly nice agreements between 
SHAM and the observation, as expected from previous work s 
dConrov. Wechsler. & Kravtsovll2006r . iTrujillo-Gomez et alfcOllh . 
We found that SHAM overpredicts the amplitude for the brightest 
sample by about 10% at all scales. This may be because we do not 
include scatter in the Viiax — M r relation. In reality, scatter be- 
tween the two quantities is expected. It has been pointed out that 
introducing scatter reduces the clustering amplitude, in particular 
for the brighter sample effectively, and makes the agreement better 
dWetzel & ; Whitdl201Ct ITrujillo-Gomez et al.ll2oTlt iReddick et all 
120121) . 

We show the HODs obtained for each binned sample in the 
bottom panel and co mpare them with th e best-fit HOD modeling 
results estimated by IZehavi et al. I J201 11). We consider the galax- 
ies assigned to the central and satellite subhalos as the central and 
satellite galaxies, respectively. The figure shows that the HODs 
of satellite galaxies of SHAM agree with the best-fit HOD model 
well. The figure also shows that the HODs of central galaxies from 
SHAM are narrower than the HOD modeling for the bright and in- 
termediate samples. It should be noted that both HOD modeling 
and SHAM give very similar predictions for galaxy clustering. 



3 METHOD FOR COLOR ASSIGNMENT 

Here we describe our method to assign galaxy color to sub- 
halos. SHAM uses the abundance as a function of subhalo 
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Figure 1. Top: The projected correlation functions for the magnitude- 
binned samples. The SHAM predictions are shown by solid lines with error 
bars ta ken from the jack-k nife resampling of eight subvolumes. The SDSS 
results l lZehavi etal.1120111) are shown by open circles with error bars. The 
amplitudes of the brightest and the faintest sample are shifted upward and 
downward by 0.1 dex for clarity. Bottom: HODs for each binned sample 
obtained by SHA M (the thick lines). For comparison, the best-fit models 
of HOD modeling fcehavi et al]|201 lb are shown by pale color thin lines. 
The solid and the dashed lines are HODs of central and satellite galaxies. 
Note that we quote the FoF mass as the distinct halo mass while the HOD 
modeling uses the virialized mass. Hence our HODs should be shifted to 
the left hand side slightly to be compared with the HOD fitting results. 



and galaxy properties to match their number densities. SHAM 
works well if the properties, e.g., galaxy luminosity and the 
subhal o maximum circular velocity, highly correlate with each 



subhal o maximum circular velocity, highly cor relate wit h eact i 
other dConrov Wechsler, & Kravtsovl l2006t iNeistein et aL1l201ll; 
lReddicketal.ll2012h . The spirit of our method is also based on 
this insight. We separate each luminosity-binned subhalo catalog 
into two groups by a secondary subhalo property. The secondary 
property is expected to be correlated well with galaxy color. Then 
the two groups correspond to the red and blue galaxies within 
the luminosity bin. In a fashion similar to SHAM, the number 
density ratio of the two subhalo groups is matched with the ob- 
served red/blue galaxy abundance ratio. We study two models 
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for the secon dary subhalo property: one is motivated by assem 
bly bias (e.gjGao. Springe & ' Whitej|2005t IWechsler et al.l | 2006 



Croton. Gap. & Wmtdl2007l : lGao & Whitj2007l : lReed et all2007 : 
Lacerna & Padillall201 If) while the other incorporates environment 
in and around subhalos. We discuss our two models below. 

3.1 Models for a proxy of galaxy color 

3.1.1 Subhalo age 

It is naive but natural to consider that the galaxy color is a proxy of 
the galaxy age and to expect that subhalo age is highly correlated 
with the age of the galax}0. We use the redshift evolution of the 
maximum circular velocity to define the formation epoch Zf olm of 
subhalos via 



1000 



Vniax (-2 — -2form) — f X 



with < / < 1 



(3) 



' max ' 

and Knax f° r the central subhalos at 2 = or 0.1, depending on the 
simulation box. Note that Zf orm of the satellite subhalos is defined 
at an earlier epoch than the accretion epoch. In practice, we identify 
the two snapshots between which / x V max is located and then we 
interpolate linearly between them to get Zf OIln . For low mass sub- 
halos or small values of /, we cannot follow down to / x Vmax due 
to the lack of mass resolution. In such cases, we define the redshift 
at which the subhalo is identified for the first time as ZformU We 
refer to this model as the age model. In this model, the subhalos 
with higher Zf orm correspond to redder galaxies. It has been shown 
extensively in the literature that th e older distinct halos are more 
clustered than the younger ones dGao. Springel. & White! 120051 : 



iGao & Whi5l2007l : lReed et al]2007l ; lLacerna & Padillal201lh" 

The choice of value of / is not unique. We find that the pre- 
dicted correlation function depends on the value of /. We show the 
impact of / in the age model in Figure [2] For example, we show 
models with / = 0.9, 0.8 and 0.6 for the -21 < M r < -20 
sample. For the fainter and brighter samples, very similar trends 
are found. The red and the blue lines represent w p or HOD for the 
subhalos with higher- and lower- Zf orm , respectively. As we expect, 
a "split" like the red/blue galaxy separation is seen in the corre- 
lation function (top panel). Furthermore, larger value of / gives 
larger split. This can be understood as the larger / value captures 
subhalos which formed at earlier epoch. Hence the predicted clus- 
tering for red galaxies is enhanced. In the bottom panel, HODs are 
shown. The peaks at M ha io ~ 2 x W 12 h~ 1 M are HODs of cen- 
tral galaxies and the power law-like curves at high halo mass regime 
are HODs of satellite galaxies. It can be seen that the satellite HOD 
is more sensitive to / than the central HOD. 

3.1.2 Local DM density around subhalos 

As our second model, we adopt the loc al DM density around sub- 
halos as the secondary subhalo property. iMandelbaum et ai1(l2006t) 
showed that early type galaxies have higher mass density pro- 
file than late type ones at 2 = ~ 0.1 in several magnitude 
bins via galaxy-galaxy lensing techniques. In particular, this mor- 
phology dependence of mass density profile is more apparent at 

3 For simplicity, we do not consider bias in age due to starburst triggered 
by mergers. 

4 As we will see later, high values of / ~ 0.9 are preferred to fit the 
observed clustering. For such /, the formation epochs of almost all subhalos 
are identified later than the first identified epoch even for the faint sample. 
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Figure 2. Impact of the adopted / value (see Equation[3j on the clustering 
and HOD in the age model. The red and the blue lines show w p or HOD for 
the subhalo subsamples with higher- and lower- Zf orm . The adopted values 
of / are 0.9, 0.8 and 0.6. Top: Projected correlation functions from the age 
model. Bottom: HODs of central and satellite galaxies as a function of host 
halo mass. 



r v > 100 /i _1 kpc. This suggests that redder galaxies tend to be 
hosted by subhalos with higher density envelope than bluer ones. 

We take a sphere with radius of Rdm centered on the center 
of a subhalo to measure the local DM density. In this model, the 
subhalos with higher local density corresponds to the redder galaxy 
populations. We take the sphere radius of ~ O(100 h~ kpc) to 
reflect our motivation from the lensing study. It should be noted 
that this property is less affected by simulation resolution than the 
age model. We refer to this model as the local density model. 

Similar to the age model, we find that the predicted correlation 
function depends on the size of sphere .Rdm- Figure [3] illustrates 
the impact of Rom in the local density model. The adopted sphere 
sizes in the figure are Rum = 500, 300 and 200 /i _1 kpc. Again, 
the predictions for the —21 < M r < —20 sample are shown. The 
red and the blue lines are w p or HOD for the subhalos with higher 
and lower local DM density, respectively. The top panel shows that 
clustering amplitudes for the subhalos with higher local density are 
higher than those for the subhalos with lower density as we expect. 
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Figure 3. Same as Figure [5] but for the local density model. The 
red and the blue lines show w p or HOD for the subhalo subsamples 
with higher and lower local density. The adopted sizes of Rom are 
500, 300 and 200 h^kpc, 



The panel also shows that the color split is larger for larger-i?DM- 
These are because that counting simulation particles in a sphere 
with radius of > 100 h~ 1 kpc leads subhalos in massive halos (i.e., 
cluster or groups sized halos) to be ones with higher local density. 
This trend is more effective for larger- Rom as shown in the HOD 
(bottom panel). The bottom panel shows that the red central galax- 
ies are assigned to higher mass distinct halos than the blue central 
galaxies. This is because that our local density measure strongly 
traces subhalo mass or host distinct halo mass itself for central sub- 
halos. 



4 RESULTS 

4.1 Color dependence of projected correlation functions 

4.1.1 Comparison with obserx'ations 

We compare our model predictions for the color depend e nce o f 
clustering with the observational results by IZehavi et al. I d201ll) . 



bright intermediate faint 

red 1.23 2.44 2.09 

blue 0.30 0.78 0.43 



Table 1. The average values of Zf orm as a function of magnitude and color 
in the age model. 

They separate the red and blue galaxies in the (g — r) — M r color- 
magnitude space with the division 

(g - r) = 0.21 - 0.03M r . (4) 

For each magnitude bin, we search the model with the lowest value 
of X 2 = Xrcd + Xbiuo. where x r 2 cd and Xbiuc are chi-square values 
for fits to the red and the blue galaxy correlation functions. 

In the age model, we take eight values of / = 
0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9 and 0.95. For the interme- 
diate sample, -21 < M r < -20, we find that / = 0.9 gives 
the least \ 2 ■ We find that, for the fainter and brighter samples, the 
model with / = 0.95 is the best model. The best model results for 
each sample are shown by thick solid lines in Figuref4] The red and 
the blue lines with error bars represent the predicted w p for the red 
and the blue galaxies, respectively. The error bars are taken from 
jack-knife subsampling with eight subvolumes. Dat a points with 
error b ars show the observational results taken from IZehavi et al.l 

boiih . 

As shown in Figure [4] our age model reproduces the observed 
color-dependent clusterings very well. The values are relatively 
high but not surprising since the time evolution of the maximum 
circular velocity is slower than mass evolution. For the Milky-Way 
sized distin ct halo, Zf orm with / = 0.95 corresponds to z ~ 1 — 2 
on average dBovlan-Kolchin et alj20ich . In TablefTJ we list the av- 
erage value of Zf OI . m in the best model for each magnitude bin and 
color. 

Next we present the results from the local density model. We 
take seven sphere sizes of i?DM = 100, 150, 200, 250, 300, 400 
and 500 h~ 1 kpc. We compare the best model shown by the thick 
lines with the observation in Figure [5] For the fainter and interme- 
diate samples, we find that the model with Rom = 200 /i _1 kpc is 
the best model. Only for the brightest sample, the best fit is given 
by the model with Rom = 250 h~ kpc. 

As well as the age model, our local density model shows 
good agreement with the observational results. It should be noted 
that the local density model is motivated by the observed lens- 
ing profiles. However the local density model systematically does 
not have enough clustering amplitude at the smallest scale, r p < 
20 ft -1 kpc, for blue galaxies. It is because that the local den- 
sity model populates less massive subhalos with more blue central 
galaxies than the age model does. We will discuss this point in more 
details below. 



4.1.2 HODs from our models 

The HODs for each magnitude bin obtained from the best mod- 
els in our two models are shown in Figure [6] The results from the 
age model and the local density model are shown by the solid and 
dashed lines. The thick red and thin blue lines represent the red and 
blue galaxy HODs, respectively. Note that the quoted distinct halo 
mass in the horizontal axis is the FoF mass. It can be seen that the 
shapes and amplitudes of HODs of satellite galaxies are very sim- 
ilar among our two models for red and blue galaxies. The primary 
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Figure 4. Comparison of the projected correlation functions between the 
best models in the a ge model (solid lines) and the observation (data points) 
IZehavi et alj201 ll) for the three magnitude binned samples. 



Figure 5. Same as Figure |4]but for the local density model. 
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Figure 6. The obtained HOD from our two models for red and blue galaxies 
in the three magnitude bins. The quoted distinct halo mass is the FoF mass. 
The results from the age model and the local density model are shown by 
the solid and dashed lines, respectively. The thick red and thin blue lines 
represent the red and blue galaxy HODs. 





bright 


intermediate 


faint 


age model 








red 


0.31 


0.41 


0.46 


blue 


0.08 


0.21 


0.23 


local density model 








red 


0.26 


0.42 


0.59 


blue 


0.16 


0.20 


0.13 



Table 2. The obtained satellite fractions from our two model for each mag- 
nitude bin. 



difference is found in HODs of central galaxies. By construction, 
the host halo mass range of central is same in the two models. In 
the local density model, however, central blue galaxies are hosted 
by low mass distinct halos which have few satellite ones. 

We here compare our two models. The most apparent differ- 
ence between them is the clustering amplitude of blue galaxies at 
the smallest scale. Systematic suppression of clustering amplitude 
can be found in the local density model. At such scale, clustering is 
dominated by signals from central-satellite pairs. As shown in Fig- 
ure [6] the local density model assigns the central blue galaxies to 
lower mass halos than the age model does whereas satellite galax- 
ies live in more massive halos. It means that the number of distinct 
halos with a central blue galaxy and more than one satellite blue 
satellites is decreased. This suppresses the clustering amplitude at 
the small scale. 

The satellite galaxy fractions for each bin from our models 
are listed in Table|2] It can be seen that the satellite fractions of red 
gala xies are higher than those of blue ones in all magnitude bins. 

IZehavi et al.1 feoill) calculated correlation functions analyti- 
cally with the model HODs. Their assumed HODs of red and blue 
galaxies have the same form but different amplitudes. They also ob- 
tained very good agreement for both red and blue samples. How- 
ever our two models do not agree exactly with their HOD forms 
though both reproduce the observed correlation functions well. 



4.2 Color dependence of lensing profile 

We have shown that both the age and the local density models can 
be tuned to reproduce the observed color dependence of projected 
correlation functions reasonably well. To discuss which model is 
better to assign galaxy color, we here consider the galaxy-galaxy 
lensing measurement AE as 



AE(r p ) = 7t (r p )£ cr = E(< r p ) - £(r p ), 



(5) 



where E is the surface mass density, y t is the tangential shear and 
S cr = (c 2 /<iivG)(D B /Di a Di) is the critical surface density (see 
iBartelmann & SchneideiEoOlL for review). We project all simula- 
tion particles around each subhalo along a dire ction in the simula- 
tion box to measure the surface mass density dMandelbaum et al 



20051: lHavas hi & White 2008; M asaki. Fukugita. & Yoshida!l201_ 
iNeistein & Khochfad^Om Using the red/blue magnitude-binned 
subhalo catalogs obtained from our models, we calculate the mean 
profile for each catalog. 

We compare the model predictions for red and blue galax- 
ies with measurements for early and late type galaxies by 
iMandelbaum et alj < l200d) . They used flux-limited samples and 
the average value of frac.deV in g, r and i bands to clas- 
sify the early and late type galaxies 0.5 for early types 
and < 0.5 for late types). The parameter frac_deV specifies the 
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galaxy flux profile as frac.deV x (de Vaucouleurs profile)+(l — 
frac_deV) x (exponential profile). As the zeroth order approxima- 
tion, we regard e arly (late) type galaxie s as red (blue) ones. Using 
the SDSS DR7 dAbazaiian et aill2009l) . we checked that 74% of 
early (89% of late) type galaxies actually lie in the red (blue) se- 
quence for —21 < M r < —20. For the fainter (brighter) sample, 
these rates are 78 (74) % and 85 (91) % for early-red and late-blue 
correspondences. Therefore the early (late) type samples primarily 
consist of red (blue) galaxies. Again , the re d/blue galaxies are sepa- 
rated via Equation l4l Tsheldon et alJ J2004 showed that early (late) 
types and red (blue) galaxies have very similar lensing profiles but 
with the different color division and early/late type definition. 

In Figure [JJ we show the results from our age model by thick 
lines with error bars for the three luminosity binned samples. The 
red and blue lines are results for red and blue galaxies, respec- 
tively. The error bars are taken from twenty seven subsamples of 
subhalos. The magenta and the cyan data points with error bars 
are the lensing prof i les of early and late type galaxies measured by 
iMandelbaum et al.1 | |2006|) . For the sample of -22 < M r < -21, 
we simply show the number weighted averages of the observed 
lensing profiles of -22 < M r < -21.5 and -21.5 < M r < -21 
samples. The used subhalo catalogs are from models with / = 0.9 
for the intermediate and / = 0.95 for the fainter and the brighter 
samples. It can be seen that the age model predictions for the red 
galaxies agree with the observational results of the early type galax- 
ies well in the three magnitude bins. However the age model does 
not trace the observation of the late type galaxies of the intermedi- 
ate and the brightest samples at 100 /i _1 kpc < r v < 1 /i _1 Mpc 
at all. 

The results from our local density model are shown in Fig- 
ure [8] in the same way as Figure [7] For the intermediate and the 
faintest bins, the adopted size for local DM density measuring is 
JioM = 200 /i -1 kpc. For the brightest bin, -Rdm = 250 /i -1 kpc 
is adopted. As well as the age model, agreements between the lo- 
cal density model and the observation are good. The most striking 
difference from the age model predictions is larger color split at 
r p 300 h~ kpc. This is because the blue galaxies in the local 
density model live in the less dense region than red ones do. The 
split is large enough to trace the late type lensing profiles. 

We discuss the difference between our two models more quan- 
titatively here. We estimate the effective mass of host distinct halos 
(M e g) and the average mass of host subhalos (M acc ) for each mag- 
nitude bin and color. (M e g) is the HOD weighted average mass of 
host distinct halos A/haio and is calculated as 



(M eff ) 



dM h 



dn (iVgal (A/halo)) „. 

>T»7 JW halo , (.0) 

aiWhalo Fi-gal 



where dn/dMh a io is the distinct halo mass function. The distinct 
halo mass A/h a io comes from the FoF algorithm. The subhalo mass 
is determined by the SubFind algorithm at the accretion epoch. We 
show the masses as a function of magnitude in Figure [9] for our 
two models. The quoted magnitude is the average value for each 
magnitude bin which can be found in lZehavi et al.1 ( 1201 ll) . The top 
and bottom panels show the results from the age and the local den- 
sity models, respectively. We see that (M e g) of red galaxies is al- 
ways greater than that of blue galaxies in both the age and the local 
density models. This is consistent with the fact that the model pre- 
dicted lensing profiles of red galaxies have higher amplitude than 
those of blue galaxies at all scales. Also in the local density model, 
(Afacc) of red galaxies is greater than that of blue galaxies. How- 
ever in the age model, (Af aC c) of blue galaxies is more massive than 
that of red galaxies. This trend can be found in HODs of central 
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Figure 7. Surface mass density profile around galaxies. The red and blue 
thick lines with error bars represent the predictions for the red and blue 
galaxies in the age model. The magenta and cy an data ponits with error 
bars are the observational results from the SDSS (Mandelbau m et alj2006h 
and show the early and late type galaxies' mass profile, respectively. The 
thin lines show results for the simulated early/late type samples (see text). 
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Figure 8. Same as Figure |7]but for the local density model. 



Figure 9. The effective mass of distinct host halos ( M c g ) and the average 
mass of host subhalos (M aC c) as a function of magnitude. {Af c g ) for red 
and blue galaxies are shown by solid and dotted lines, respectively. (M aC c) 
for red and blue galaxies are shown by dashed and dotted-dashed lines, 
respectively. Top: The results for the age model. Bottom: Same as Top panel 
but for the local density model. 



galaxies from the age model (see Figure [6j. This is because more 
massive halos have more rapid mass growth on average (see e.g., 
IWechsler et al .l2002h . Hence, with our definition of the subhalo for- 
mation epoch, more massive subhalos tend to be blue galaxies in 
our age model. Thus the resulting low (high) mass of host subha- 
los of red (blue) galaxies in the age model makes the color split of 
the lensing profiles smaller than in the local density model even at 
r p < a few x 100 /i _1 kpc. The characteristic masses are listed in 
Tabled 

To discus s further, we make a more " fair" comparison with the 
observation of iMandelbaum et alj d2006l) . We match the red/blue 
fraction of our simulated galaxies in a magnitude bin to the fraction 
we obtain from the actual SDSS galaxy catalog for the early/late 
type sample. To this end, we blend the red/blue galaxies randomly. 

The results of lensing measurements from the matched 
early/late type galaxy catalogs are shown by thin lines in Figures 
|7] and [8] without errorbars. The red and blue thin lines represent 
lensing profiles of the early and late types respectively. We see the 
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Table 3. The effective average mass of host distinct halos (M 8 g) and the 
average mass of host subhalos (Afacc) from our models for each magnitude 
bin and color. The unit is /i Mq. 

effect of the color fraction matching clearly at r p > 300 /i _1 kpc. 
Figure[8]shows fairly good agreement between the predictions from 
the modified subhalo samples with the local density model and the 
observations. However the agreement is less impressive for the age 
model subhalo samples, as can be seen in Figure [7] This implies 
that the local DM density is more crucial for the galaxy color as- 
signment model than the subhalo formation epoch. 



5 SUMMARY AND DISCUSSION 

We have extended the subhalo abundance matching (SHAM) to de- 
velop a new scheme that assigns galaxy color as well as luminosity 
to subhalos. We consider the subhalo age and the local density as 
a secondary subhalo property which is expected to be correlated 
with galaxy color. Technically, we divide a magnitude-binned sub- 
halo catalog obtained by SHAM into two samples by the secondary 
property. The two samples are then meant to represent red galax- 
ies and blue ones. In a similar fashion to SHAM, the abundance 
ratio of the red and blue galaxies is matched to the observed ratio. 
We have studied the spatial clustering of the red and blue galaxies 
based on the subhalo age model and the local density model. Over- 
all, the two models re produce the observ ed color-dependent galaxy 
clustering properties dZehavi et al .11201 ll) reasonably well. It is en- 
couraging that SHAM can be extended successfully in this way by 
using a secondary subhalo property that can be easily measured. 

The clustering amplitudes of the blue galaxies in our local 
density model is systematically smaller than the observation at 
r p < 20 h~ kpc, where the central-satellite signal dominates. 
This is because, in the local density model, central blue galaxies 
tend to be hosted by relatively small distinct halos (see Figure [6). 
The number of pairs of a blue central and a blue satellite is then 
smaller than in the observation. The small-scale clustering can be 
better matched to the observation by introducing scatter in the lo- 
cal DM density-color relation. In the real universe, there are sub- 
stantial scatters in the relation. It is straightforward to model the 
scatters by perturbing the measured local DM density with a cer- 
tain probability distribution, e.g., a Gaussian with a finite variance. 
Such an 'improved' model would then assign a blue central galaxy 
more likely to a massive distinct halo. We have explicitly checked 
that introducing scatters in the local density - color relation, in the 
above ad hoc manner, improves the agreement in the small-scale 
clustering of the blue galaxies. 

We have examined a few other subhalo properties as an alter- 
native proxy for galaxy color, for instance, the density concentra- 



tion and the spin parameter. We have also tried with slightly differ- 
ent definitions for the concentration, spin, age and local DM den- 
sity. None of them yields better agreement with the observed cor- 
relation functions. We thus conclude that our original age model 
and the local density model are the best among those we have ex- 
amined. It is worth discussing further the success of our local den- 
sity model. One can naively expect that galaxies in higher density 
environment formed earlier and that they are not forming stars ac- 
tively any more. Thus the local matter density could be related to 
the galaxy color. The local DM density around a subhalo within 
100 /i _1 kpc-order sphere closely traces the subhalo mass or, for 
central galaxies, the mass of the host parent halo. Then the sub- 
halo mass itself could be used as a proxy for the galaxy color rather 
than the local density. We can test the idea by using the subhalo 
masses determined by the SubFind algorithm at the observation 
time. In this model, more (less) massive subhalos correspond sim- 
ply to redder (bluer) galaxies similarly to the local density model. 
We have run the same calculations in the previous sections. The 
subhalo mass model does not reproduce at all the observed color- 
dependent clustering. Because subhalos in a massive distinct halo 
typically have small masses, owing to tidal mass stripping, they 
tend to be assigned a blue galaxy in this model. Indeed, the se- 
lected blue galaxies show strong clustering at all length scales. We 
argue that some environmental effect, in terms of local density, is 
important for color assignment for the satellite galaxies. 

It is crucial to test our model predictions against not only spa- 
tial clustering but also other measurement(s). We propose to use 
the color dependence of l ensing profile, motivated by the observed 
morphology dependence ( IMandelbaum et al.l200d) . To this end, we 
regard the early (late) types as the red (blue) galaxies and calculate 
the lensing profiles for the two populations. The predictions from 
our two models agree with the observed red galaxy lensing profile 
reasonably well. In particular the local density model reproduces 
the observed lensing progile for the blue (late type) galaxies. How- 
ever, while the local density model shows a large color (morphol- 
ogy) split at 100 h~ kpc < r p < 1 /i _1 Mpc, the age model does 
not account for the low amplitude of the lensing profiles of the late 
(blue) galaxies with -21 < M r < -20 and -22 < M r < -21 
at 100 /i _1 kpc < r p < 1 /i _1 Mpc at all. 

The fact that the age model does not reproduce the observed 
color (morphology) dependence of lensing profiles might indi- 
cate that galaxy color and subhalo age are not tightly correlated 
with each other. There is an ob servational hint for this notion. 
iTinker. Wetzel. & Conrovl d201 ll) studied the fraction of central 
galaxies that are not forming stars actively ("quenched centrals") as 
a function of the local galaxy density using a SDSS group catalog. 
They show that the fraction of passive galaxies does not agree with 
the fraction of old halos, in particular for relatively small galax- 
ies with log[M„ / (h~ 2 Mq)] < 10. Their findings are qualitatively 
consistent with our result. There is not a tight relation between the 
galaxy color and the age of the host (sub)halo, or one may need to 
devise a rather complex conversion from one to the other. We can 
think of using the subhalo assembly history of red and blue galax- 
ies directly obtained from A^-body simulations. Further studies are 
clearly needed. 

The abundance matching is a powerful technique to populate 
a cosmological simulation with galaxies. It is encouraging that our 
simple models reproduce the observed color-dependent clustering 
well. There are a variety of choices of the secondary subhalo prop- 
erty as a proxy of galaxy color, and thus it is important to test mod- 
els using other observations independent of spatial clustering, as 
we have done in the present paper. 
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